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1. Introduction 



Considerable tiieoretical and experimental attention has been devoted to the extrac- 
tion of the Cabibbo-Kobayashi-Maskawa element \Vcb\ from both exclusive and inclusive 
semileptonic decays of the B meson. Because of the key role this parameter plays in the in- 
vestigation of rare decays and CP violation in the third quark generation, it is important to 
make as many independent determinations of \Vcb\ as possible. The baryonic semileptonic 
decay A;, A^^z^ provides an opportunity to extract \Vcb\ in a fashion as theoretically 
clean as in S ^ Dii? or B D*iiJ, because heavy quark symmetry predicts for both 
cases a single universal form factor in the heavy quark limit 0] , including the normalization 
of 1 + 0{l/ml) at zero recoil||,|. 

The prospect for experimentally determining the dependence of A;, —>■ A^^z/ form 
factors is promising, because the method applied by ALEPH[0 to B D*iu should work 
equally wellg for baryons. The cumulative data sets of ALEPH, DELPHI, and OPALQ 
contain about 250 A^ semileptonic events; an additional 200 events have been observed at 
CDF[|TD[. This raises the possibility of extracting |Vc&| from baryon decays. 

To fully exploit this possibility, a model-independent parametrization of the A;, 
Aciiy form factors is desirable. This is important because the kinematic rate vanishes at 
zero recoil, so an extrapolation of the data based on some parametrization is required. 

In this paper we apply the parametrization imposed by QCD on the physical form 
factors. The requirement of compatibility with QCD is obtained through the application of 
dispersion relations based on the analyticity properties of form factors as functions of the 
momentum transfer variable. The method of extracting information on amplitudes by this 



form of complex analysis is quite old[|TT|, and has been applied to the study of semileptonic 
decays of light mesons in a more contemporary language in Ref. [0. Its application to 
heavy quark systems has received much attention in more recent years |[T^- [T9|] . 

In general, these applications say little about the overall normalization of form factors, 
but encode a great deal of information about the shape. The expression of this information 
manifests itself in a simple parametrization[T^,T^ that spans the functional space of form 
factors allowed by QCD dispersion relations. In this work we apply the parametrization, 
supplemented by a new development presented here, to the decay A^ Acify. 

The technical development introduced for A;, Aciu also has important implications 
for form factors in meson decays such as B ^ D^*^iiy and B niiy. It allows one to obtain 
form factor parametrizations with smaller uncertainties and fewer parameters than those 
discussed in earlier works [|1^,|T^ . 
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The paper is organized as follows: In Sec. 2 we outline the derivation of form fac- 
tor bounds from QCD dispersion relations and introduce the technical improvement that 
strengthens our constraints. We construct the form factor parametrization that obeys 
these bounds, and discuss the inclusion of heavy resonances, in Sec. 3. In Sec. 4 we dis- 
cuss the numerics of the heavy resonance masses and obtain information on the quality of 
parametrizations for Aj, Aciu form factors. Section 5 describes the implications of our 
technical improvement to meson decays. We conclude in Sec. 6. 



2. Dispersion Relation Inequalities 

The QCD matrix elements governing the semileptonic decay A;, - 
expressed in terms of form factors defined by 



Ac^z^ may be 



(2.1) 



where v = p/Mi>^^^ and v' = p' /M^^ are meson velocities, = c'j'^b and A'^ = c'j'^'j^b are 
polar and axial vector flavor- changing currents, and the form factors are functions of the 
momentum transfer q'^ = {p — p')"^- In terms of these form factors, the differential decay 
width for A^ Kci-i' with a massless charged lepton I is 



.{(Q_) [2q2|^^|2 ^ M^lif^p] + (Q^) [2q2|^^|2 ^ M^lif^p] }, 

(2.2) 
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where the terms proportional to and alone give the partial widths for trans- 



versely polarized intermediate W bosons [POI , whereas 
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(2.3) 



are the partial wave amplitudes for a longitudinally polarized VF's, with 



(2.4) 



The factor 



(2.5) 



is an invariant kinematic quantity related to the three-momentum of the in the 
rest frame of the decaying Af,, and M,m = M\^, M\^, respectively. The form factor 
combinations — ^) and — which appear with the Lorentz structure q'^, give 
contributions to the rate proportional to the lepton mass squared mj , and are consequently 
mainly of interest for constraining models. 

We begin our derivation of constraints from dispersion relations following the well- 
known methods developed by the authors listed in Refs. [|TT| , p!^ . In QCD, the two-point 
function of a flavor-changing current J = V,A, or V — A, 

W/iq) = {qf^q" - q^ g'"')U^ (q^) + g'^'U^iq^) = i j rfic e*'?^(0|TJ'^(a;) J^^(O) |0), (2.6) 

can be rendered finite by making one subtraction, leading to the dispersion relations 

TL. 2n dlfi'^ 1 , Imn'^'^(t) , , 

The functions Xj'^ici^) may be computed reliably in perturbative QCD for values of q^ 
far from the kinematic region where the current J can create resonances: specifically, 
{nib + Wc)Aqcd ^ {rrib + TTicY ~ 1^ ■ For resonances containing a heavy quark, it is 
sufficient to take q^ = 0. At one loop, Xv ~ '^■'^ ' 10""^ and Xa ~ 2-2 ■ 10~^. We also make 
use of the quantity Xj = Xj (0) ~ I af^'XjlO), corresponding to the combination of 11 j and 
Hj that gives Hj at q^ = 0. The full one-loop expressions for q"^ — are 

Xv(^) = X^(-w) 



327r2m2(l -^2)5 

X [(1 - u^){3 + 4u- 21u^ + 40u^ - 2\u^ + + 3w^) 
+ 12^2(2 -3w + 2^2) inw^], 

Xv{n) = x\{-n) = g^2(/_^2)3 [(1 + « + " + v?) - GtiMnw^] , 

where u = rric/rnb is the ratio of quark masses. For bottom and charm quark masses of 
nil, = 4.5 GeV and nic = 1.5 GeV, the one-loop values of M'^xv and M'^xa are 1.5 ■ 10~^ 
and 9.0- 10"'^, respectively. One should keep in mind that 0{as) corrections to the bounds 
derived below enter as corrections to these values. 

The absorptive part Imnj'^(q'2) is obtained by inserting on-shell states between the 
two currents on the right-hand side of Eq. ( |2.6| ). For any four- vector n^, the quantity 
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*Tl^j^ni, is a sum of positive-definite terms, so one can obtain strict inequalities by con- 
centrating on tlie term witli the two-particle intermediate state A^Ac- The contribution of 
A^Ac pairs to the right-hand side of ( P77| ) enters as 



n*Imn^^_^n^ > 9{q^ - {M + mf)^-^-= j dQp ^MmXi{n* ■ n) 



m 
M 



—X^\{q-p) ■ 
m 



(2.9) 



-f X5{p ■ n*){q -p) ■ n + X^{p ■ n){q- p) ■ ?7,*|, 



where 



Xi = (u;-l)|Fi|2 + (l + w)|G'i|2, 

^2 = (1 + oo) |F2|2 + {F*F2 + FiF;) + {uj-l) 1^21' + {GIG2 + GiGl), 
Xs = {F*Gi + FiGl), 

X4 = (1 + u) \Fs\^ + iF*Fs + F^F;) + (u - 1) \Gs\^ - (G^Gs + GiG*), 

X5 = (1 -f u) f;Fs + - 1) G*2Gs + \Fi\^ + FiF* + F*Fs + \Gi\^ - G^Gl + G^Gs, 

(2.10) 

with uj = (M^ + — q"^) / {2Mm), and the integration is over all directions of the mo- 
mentum vector p. For massless leptons the differential width Eq. ( p.2|) is proportional 
to the space-space components Imllj of the two-point function, so the choice n = (0, n) 
leads to inequalities on the physically interesting form factors Fi, Gi, Hy, and Ha- For 
completeness, we also consider n = (1,0), which leads to inequalities on the combinations 



Fo = — 



(M - m)Fi + -i-(g2 ^ ^2 _ ^2^^^ _J_^q2_M^ + m')F^ 
ZM 2m 

(M + m)Gi - -i-(g2 ^ m2 - m^)G2 + ^^{q^ - + m^)G3 
2M 2m 



(2.11) 



However, since can be an arbitrary g^-dependent four-vector, there is considerable 
freedom to constrain other combinations of form factors. 

The analysis simplifies when Eq. (|2.7| ) is written in terms of the conformally trans- 
formed variable z defined by 



l + z {M + my-q^ 



(2.12) 



l-z V AN Mm V 2N 

This expression differs from that previously used in the literature by the inclusion of 
the factor A^; the two agree when = 1. Upon choosing the principal branch of the 
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square root in this expression, the change of variables ^ z maps the two sides of the cut 
> (M+m)^ to the unit circle \z\ = 1, with the rest of the plane mapped to the interior 
of the unit circle (Fig. 1). In particular, the real values — oo < q^ < {M—m)'^—4{N—l)Mm 



and (M 



4(A^ — l)Mm < q^ < (M + m)^ are mapped to the real axis, 1 > z > 



and > z > —1 respectively. Physically, this means that the kinematic region relevant to 
the process vacuum Mm, q^ > [M + m)^, lies on the unit circle, while the region for 
semileptonic decay M miu for massless i, < q^ < (M — m)^, lies inside the circle on 
the real z-axis. Specifically, 



(l-^)2_2^(v/iV-l) 
(l + v^)2 + 2v/f(ViV-l) 



^2.13) 




Figure 1. Important kinematic points in the z-plane for the related processes (vacuum 
Mm) and M miu using the map Eq. 0^.1^ ) ■ a) (z = —1): The branch point 
q^ = {M + m)^, threshold for (vacuum Mm). The two sides of the branch cut 
q^ > (M + m)'^ map to the upper and lower halves of the unit circle, with q^ — > 
-\-oo corresponding to z — > +1 along the circle, b) Location of resonant poles below 
the (vacuum Mm) threshold. c,d) Zmim Zmax defined in Eq. (2.13). The segment 
(Zmim Zmax) is the kinematic region for M — > mii^, and N is chosen in Eq. ( 2.12 ) so 
that z = lies in this interval, (f' ^ —oo corresponds to z — > +1 along the real axis. 

where r = m/M. Inasmuch as both l^maxl, kmini ^ 1, semileptonic decay possesses a 
small kinematic expansion parameter. One can use the free parameter to adjust where 
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on the real axis the interval (^mim ^max) lies, but it is clear that allowed z values are 
uniformly smallest when z^^^ < < ^max, for which 

1 < < (M + m)V(4Mm). (2.14) 

Choosing saturates the last remaining degree of freedom allowed by the Schwarz- 
Christoffel transformation between the cut plane and the unit disc. may be kinematically 
interpreted as the value of (1 +u;)/2 for which one obtains z = in Eq. ( |2.12| ). We explore 
the consequences of varying in the next section. 

Written in terms of z^ the inequalities from Eqs. (|2.7| )- (|2.10|) now read 



(2.15) 



The contour C is the unit circle. The weighting functions (^i are constructed to be analytic 
inside the unit circle by multiplying kinematic factors by functions that are phases on the 
unit circle. For example, 



M2 
M2 

M2 



[(1 + r)(l - + 2ViV^(l + z)fl{\ - zf, 
4r 



1 -Z)2 

(1 + ^)2 



[(VAT - l)z + (VAT + 1)] 



4rA^ 



(1-^)2- 

After some algebra, one obtains the weighting functions 

[(1 + r)(l - 2) + 2ViV^(l + [(ViV - \)z + (TiV + 

where p, and s depend on the form factors Fi as listed in Table 1. 



(2.16) 



(2.17) 



i 


F^ 


1/k 


P 


s 


1 

2 
3 
4 
5 
6 


Fo 
Fi 

Hv 
Go 
Gi 
Ha 


GnM^Xv 
12nM^Xv 
^t^Xa 

GnM^XA 
UnM^XA 


1 




1 
1 





1 




1 



Table 1. Factors entering Eq. ( \2.17\) for the form factors Fi. 
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3. Parametrization of Form Factors 

Equation ( 2.15 ) constrains the form factors Fi along the unit circle in the complex z 
plane, but for semileptonic decay we are interested in the physical region along the real 
axis, ^min ^ z < -Zmax- If the product (f)i{z)Fi{z) were analytic for 1^1 < 1, one could simply 
perform a Taylor expansion about z = 0. While the weighting functions (pi are analytic 
inside the unit circle, the form factors Fi are not. Contributions from intermediate states 
with masses below the A;, Ac threshold lead to cuts and poles in Fi{z) for \z\ < 1. While we 
expect the contribution from cuts to be unimportant [p!8| , the more singular nature of poles 
requires a careful treatment [jl^ . The use of Blaschke factor s |]2T[] permits one to eliminate 
pole contributions fl^ given only their masses a/^ = rrij, which are converted to positions 
Zj inside the unit circle via Eq. ( p.l2|) . The Blaschke factors are defined by 

{z~Zj) 



V 



(1 - ZjZ) 
-L.-L (1 - ZjZ) 

where V and A indicate the poles appropriate to the polar vector and axial vector form 
factors in Eq. (|2.1| ), respectively. Resonances with spin-parity = 0"*", 1~ couple to form 
factors of the polar vector current, while resonances with = 0~, 1"*" couple to the axial 
vector current. The form factors coupled to spin-0 resonances are exactly those that are 
multiplied by the Lorentz structure g^, and consequently only appear in the differential 
width multiplied by the factor mj. It follows that, in the massless lepton limit, only polar 
and axial vector resonance masses are needed in Eq. ( |3.1j ). 

The usefulness of the Blaschke method lies in the fact that each factor (z — Zj)/ (l — zjz) 
serves to eliminate the resonance pole behavior l/{z — Zj) of the jth pole, but is unimodular 
for \z\ = 1. This holds also for any product of these factors. Therefore, since each Pi is 
unimodular on the unit circle, one may replace (piFi with PiCpiFi in the bound Eq. ( |2.15| ) 
without changing the result. Since now both [PiFi] and are analytic on the unit disc, 
one can perform a Taylor expansion about ^ = 0, 

^ oo 

Substituting this expression into the modified Eq. ( p.l5|) gives 



Y,\an?<l. (3.3) 



n=0 
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It should be pointed out that one pays a price for these factors. The fact that one can 
ehminate pole contributions without reference to their residues means that this method 
must apply equally well for all allowed values of the pole residues, and therefore the bound 
of Eqs. ( p.2| ) and ( |3.3| ) is necessarily weaker than the corresponding bound obtained if 
one knew anything about the sizes of the individual residues. In fact, the bound derived 
above is not very restrictive at all, until the constraints supplied by fitting to one or more 
experimental data points are imposed. Then the range of allowed form factors becomes 
surprisingly small, as was showed in Refs. [|T7| , p!8[| . 

Another strength of the parametrization Eq. ( |3.2| ) is that, when r = m/M is not 
much less than 1, the scaled kinematic variable z defined in Eq. (|2.12| ) is quite small. For 
example, if = 1 the decay Kclv has Zmin = and ;2max = 0.049. Using Eq. ( |37^ ) 

it is clear that the series in Eq. (|3.2|) converges quickly, and one can define an uncertainty 
from truncating the series after only a few terms. It follows that one obtains a functional 
form for the form factor under consideration (Eq. ( |3.2| )) over the whole kinematic range 
for semileptonic decay, with a well-controlled and small uncertainty (because z is small); 
only the first few parameters a^, each of which lies in a restricted range (Eq. ( p.3|) ), are 
needed. 

One can do even better by allowing for the parameter A^ 7^ 1. To illustrate this, let 
us define an approximation by truncating after the Qth term: 



n=0 



Then the maximum error incurred by truncating after Q terms is just 
max|F,(z) - Ff^{z)\ = max |a„z-| 



< max 



1 



oc 



J2 (3-5) 



< max 



2€(2„iin,2max) \Pi{z) (l)i{z) \ y/ 1 - ^2 



where we have used the Schwarz inequality and Eq. (|3.3| ). Inasmuch as \Pi(z)(f)i{z)\\/l — z^ 
varies slowly over the small interval Zmin ^ z < Zmax, the truncation error is driven 
by 1^*5+^1. One then sees that the minimum error occurs for a value of A^ such that 
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z^in — —Zmax- This simple observation allows for an improvement over the corresponding 
expression in Refs. which used = 1 and z^i^ = 0. The interval (^max — z^in) is 

nearly independent of A^, so -Zmax for optimal values of is approximately half of its value 
for = 1. Although 2; is a function of A^ through Eq. (p.l2| ), the Blaschke products Pi{z) 
for fixed and pole masses are independent of A^, while the weighting functions (p{z) are 
only weakly dependent on A^. It follows that for optimal A^, truncation errors are reduced 
by a factor of nearly 2*^+^, which is appreciable even for small Q. 

The reader should be reminded that the functional form Eq. ( |3.2| ) has known, bounded, 
and small truncation errors, regardless of the specifics of the experimental data. A form- 
factor expansion in (u; — 1) truncated after linear or quadratic order introduces theoretical 
uncertainties that can be substantial fl^, and unlike ( |3.2|) is not a priori guaranteed by 
theory to fit the data. 



4. Heavy Resonances 

To complete the description of the A;, Kcli' form factors provided by Eq. (|3.2|), 
one must construct the Blaschke products Pi{z). The relevant poles are resonances below 
the AfoAc threshold that couple to the b ^ c currents. These are the polar and axial vector 
Be mesons with masses below Ma^ + ~ 7.9 GeV. Although results on Be mesons are 
still very preliminary ||22|| , their masses are reliably calculable by interpolating between the 
charmonium and bottomonium spectra. The quark model in a confining potential predicts 
a number of relevant Be polar vector states with the spectroscopic quantum numbers 
^Si and ^Di. In the bottomonium nomenclature, the analogous states are called T{nS) 
and T{nD), respectively. There are also relevant axial vector states, which divide into 
nearly degenerate pairs of ^Pi and "^Pi states broken only by spin-orbit and hyperfine 
splittings. The bottomonium analogies are respectively called hbijiP) and Xbi(j^P)- The 
values obtained by various researchers [p3| -p5| using different approximations in potential 
models agree within a few MeV. We use the values of Ref. , supplemented by values for 
a few higher states above the BD^*^ threshold but below the A^Ac threshold ||26|| . Although 
the masses of the latter states do not yet appear in the published literature, their calculation 
is no more difficult than that of those already published; however, one should keep in mind 
that the simple nonrelativistic bound-state picture for such B^ resonances is no longer fully 
justified: The coupling of open channels to the resonances above the BD threshold can 
modify the details of the mass spectrum. Nevertheless, experience from the J/ij} and T 
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systems shows that such calculations give a reasonable account of masses above the DD 
and BB thresholds. We present the B^ pole mass values used in Table 2. 



Type 


Source 


Masses (GeV) 


Vector 


Ref. 
Ref. 


24 1 
26| 


6.337, 6.899, 7.012, 7.280 
7.350, 7.594, 7.646, 7.872 7.913 


Axial 


Ref. 
Ref. 


24 1 
26| 


6.730, 6.736, 7.135, 7.142 
7.470, 7.470, 7.757, 7.757 



Table 2. Assumed Be pole masses used in this work. A complete set of masses 
for either parity consists of the values from Ref. together with the newly- 
computed values from Quiggp^]. 

The truncation errors on our parametrizations depend on the Blaschke factors only 
through the endpoint values ^(^min) and P(2max), because empirically the maximum trun- 
cation errors from Eq. ( |3.5| ) occur at either z = Zniin or Zmax- These values of P do not 
depend on A^, but they do depend on the masses of states between the BD and A{,Ac 
thresholds. We present the truncation errors for Q = 2 (a two-parameter description if 
given the overall normalization, or three parameters otherwise) in Table 3. Truncation 
errors for Q = 3 are roughly a factor of 40 times smaller because of the additional sup- 
pression of about /2 described in Sec. 3. To compare these absolute truncation errors 
to the expected size of the form factors, note that heavy quark symmetry predicts the 
normalization of Fi and Gi at zero recoil to be unity, while the other form factors are 
0{l/m,c) in this limit[Q. Empirically, the value N = 1.09 minimizes the truncation errors, 
although this number may change slightly depending on the actual values of P(2;min) and 
-P(2max)- These truncation errors are substantially smaller than those for A = 1, which is 
28% for Fi, for example. Randomly altering the resonance masses by 20-50 MeV alters the 
truncation errors presented in Table 3 by less than one part in ten in each case, indicating 
the insensitivity of the errors to the particular assumptions of the potential model used to 
compute Be masses. 





(Max. Trunc. error) ■ 10^ 


Aideal 


Fo 


4.0 


1.088 


Fi 


5.1 


1.088 


Hv 


20.0 


1.089 


Go 


4.3 


1.089 


Gi 


6.6 


1.088 


Ha 


25.7 


1.089 
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Table 3. Truncation errors from Eqs. ( ]3.4| ) and ( ]3.5| ) of the Ab — * Ac^z? form 
factors for the pole masses given in Table 2 with Q = 2. A^ideai is the value of 
N for which these errors are minimized. Including all other known corrections 
multiplies the truncation errors by the coefficient B <1A described in the text. 

Finally, we briefly commented in Sec. 3 that the effect of branch cuts due to non- 
resonant on-shell intermediate states connecting the vacuum to A^Ac are not numerically 
significant, just as in B D^*Up. In Ref. M\ their effects were accommodated by re- 
laxing the bound (|3.3| ) to E|a„p < B^, with B < 1.05. In the present case, however, the 
A^Ac threshold is higher than the BD^*^ threshold, and so many more cuts are possible. 
To argue that their total effect is still negligible requires a more detailed analysis than 
appears in 
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We begin with the observation that multiparticle intermediate states connect the 
current to A;,Ac, and thus appear as loop diagrams. Each additional intermediate state 
introduces a loop and large- A^coior factor of l/( IGtt^A^c^^), so we limit ourselves to two- 



particle intermediates. In reference [^, we modeled the contribution of such a state with a 
square-root cut, as follows. If the combined mass of the intermediate particles corresponds 
to z = Zcut, the cut contribution to B is 



5B 



2-K 



^ / de\g,^,{z)(j){z)\' 

ZTT 



1/2 



(4.1) 



where z = e*^, (t){z) is the weighting function, and 



9cv.t{z) 



)(1 



ZZcut J 



l + z 



Ma, \ (l-2)(l-2cut) 



2(1 



(4.2) 



models the cut. 

The coupling constant c may be written as c = fg/8n, where / is the coupling between 
the current and the two-particle intermediate state, g is the overlap between these states 
and A^Ac, and the Stt comes from the loop. 

Two-particle intermediate states with strangeness are OZI suppressed, so the only 
kinematically allowed possibilities are B — D resonances and states with one Be resonance 
and one light unfiavored meson. We may further limit our consideration to S-waves only, 
since higher partial amplitudes are suppressed by powers of the available three-momentum, 
O(Aqcd) because we are in the resonance region, divided by the center-of-mass energy, 
0{Mb + Md). Furthermore, two-particle intermediates in which one of the particles has 
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strong isospin (such as BcTt) are heavily suppressed by isospin-breaking coefficients, and 
will be ignored. 

With these restrictions, the decay B D*£i/ admits only two very short cuts in the 1~ 
channel {Bc{^Pi)r] and Bc{^Pi)r]), with Zcut ~ —0.88, and three in the 1"*" channel, with the 
lowest being B*rj with Zcut = —0.48. There are anomalous cuts starting at z ~ —0.35, but 
they are proportional to a B*-B-n coupling g'^ that is probably quite small|]T^. B Dii> 



has no cuts in the relevant 1 channel, anomalous or otherwise. These branch cuts lead 
to even smaller uncertainties for B D^*^li> than anticipated in Ref. On the other 



hand, once the threshold is raised to the A;, + Ac mass, many more cuts occur. We count 
32 in the 1~ channel and 28 in the 1"*" channel, including such exotic combinations as 
i?c(2 '^Po(6700)) + /ii(1170). To justify our neglect of so many cuts requires a more careful 
study. 

The great majority of cuts in either channel arise from a B^. resonance and a light 
unfiavored meson. To get a rough estimate of their combined contribution, we take them 
all to have the same coupling c, compute the contribution to B of each state separately, 
and add them in quadrature (to simulate random phases between the various cuts). We 
find that B is only increased to 1 + 0.5c for the 1~ channel and 1 + 0.6c for the 1"*" channel. 
Furthermore, we expect c = fg/Sir to be very small. To judge the typical size of /, for 
example, we can extract the coupling /^^ of a vector current to the J/ifj-rj state from 
charmonium radiative decay data. Using a phenomenological Lagrangian term 

6C = ieUrjM^ij^A^T) (4.3) 

gives a coupling /^^ < 10~^, while bottomonium data gives an even smaller bound on frrj- 
Taking the interpolated value of fs^rj as representative, we conclude that cuts involving 
Be resonances give negligible contribution to B. 

On the other hand, the intermediate states consisting of a B-D pair may have much 
larger couplings, with fg conceivably of order one, but the cuts begin much closer to the 
A^Ac threshold and are fewer in number (six for the 1~ channel, with minimum Zcut = 
—0.55, three for the 1"*" channel, with minimum Zcut = —0.39). Added in quadrature, the 
bounding parameter B in the vector (axial) channel receives a 0.12c (0.18c) correction. 



However, since our uncertainty estimates to obtain i?[|T8[ are added in quadrature, the 
conservative total correction to Eq. ( |3.3| ) given by S < 1.4 remains unchanged in this 
work, even if c = 0{1). In addition to multi-particle cuts, these uncertainty estimates 
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include corrections to the heavy quark hmit, choices of quark mass values, and 0{as) 
perturbative corrections, as well as others. A value of B = 1.4 has the effect of increasing 
all of our truncation errors by a factor of 1.4; such increases on a few-percent truncation 
error do not hamper our conclusions. 



5. Implications for Meson Decays 

The freedom to choose in Eq. ( p.l2|) has significant implications for semileptonic 
meson decays like B D'^*^Iv and B ulu. For B Dlu and B D*li>, the form 
factors defined by 



{D%p', e)\A^ \Bip)) = foe*^ + {e* ■ p)[a+ip + pT + a_(p - pY] 

{D{p')\ \B{p)) =u{p+p'r + f-ip-p'r 



(5.1) 



and 



1 

m 



2q'k'a+--iq'-M' + m')fo 



(5.2) 



where now M = Mb and m = Mj^{*) , have weighting functions 

[(1 + r){l-z)+ 2Vn^{1 + z)]-^'+P\{y/N - l)z + (ViV + 1)]* ; 



(5.3) 



where n/ is an isospin Clebsch-Gordan factor that is 2 for B — > D^*^!^ and | for B txEv, 
r = m/M, and the values of k, p, and s are given in Table 4. The parameters p and s here 
do not have precisely the same origin as those in Eq. ( |2.17| ). Equation ( |5.3|) agrees with 



Ref. 181 in the limit ^ 1. 



i 




1/k 


p 


s 





fo 




1 


3 


1 


Fi 


24nM^XA 


1 


4 


2 


9 


12nM^Xv 


3 


1 


3 


/+ 


6nM^Xv 


3 


2 



Table 4. Factors entering Eq. ( \5.3i ) for the meson form factors Fi. 

From the discussion in Sec. 3, one expects that by choosing the optimal value between 
A^ = 1.10 and 1.12 (the exact value depends upon the form factor), the truncation errors 
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for a three-parameter description should be reduced by roughly a factor of 8 over the 



= 1 values in Ref. [|T^, while the truncation errors for a two-parameter fit should be 
reduced by a factor of about 4. From direct calculation, these numbers are 6.1-7.2 and 
3.3-3.7, in accord with expectations. Since one of the parameters is determined from the 
normalization of the form factor at zero recoil by heavy quark symmetry, the ~ 1.1 
choice provides a one-parameter description of each of the form factors with truncation 
errors (relative to the normalization at zero recoil) of 6.4%, 5.0%, 2.9%, and 30.9% for 

/o, g, and Fi, respectively. For the form factor g, for example, the conclusion is 
that using the normalization at zero recoil (which determines ag) plus the slope at the 
same point (which determines ai), one obtains the shape of g over the entire kinematic 
range with a theoretical error of no more than 2.9%. Uncertainties like corrections to 
heavy quark symmetry, perturbative QCD corrections, and so on are expected to increase 
this theoretical error to no more than 4.1% and possibly much less[0. For example, the 
truncation error vanishes at zero recoil, where the normalization is known. 

We exhibit the possible shapes of g{uj)/g{l) in Fig. 2. The lowest curve corresponds 
to the saturation of Eq. ( |3.3|) by ai = —1, while the top curve corresponds to ai = 0.115, 
the largest value allowed by the the Bjorken inequality on the slope of the Isgur-Wise 
function improved by 0{as) corrections]^^. Intermediate curves correspond to equally 
spaced values of ai. 

The form factor defined by replacing the D meson in Eq. ( ^.1|) with a pion has 
the same weighting function (f)^ as f^, but with the D mass replaced with the pion 
mass and xv = 9.5 • 10~'^/m^. The Blaschke factor for the pion form factor is = 
{z—z^)/{l—zz^), with z^ corresponding to the B* mass, the sole resonance below threshold. 
One might expect a suitable choice of to lead to \zmin\j \zmax\ of about half the = 1 
value of z^"^^ ~ 0.5. Unfortunately, one cannot ignore the variation of P{z)(j){z)y/1 — z"^ 
over the kinematic range, so the reduction in truncation errors is less spectacular than in 
the B ^ D system. It is nevertheless a significant improvement: Where for A^ = 1, a 
four-parameter description implied an absolute truncation error of 1.35, for A^ = 2.1 the 
absolute truncation error is 0.37. A five-parameter description yields, for A^ = 2.27, a 0.11 
truncation error. Even for a three-parameter description, the absolute truncation error at 
A^ = 1.86 is 1.19 over the large kinematic range of B Triu. In typical models, f^{q^) 
varies from around 0.3 at = 0[|2^ to order 10 at = {Mb — m7r)^[|l6|- Other models 
may give smaller values of /+(?max) ^^^^ rapidly away from g^a,x' then it is useful 
and straightforward to choose A^ in order to minimize the truncation errors relative to 
particular points fixed by a model, instead of the absolute errors quoted above. 
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6. Conclusions 



The analysis of dispersion relations and analyticity properties of strong-interaction 
amplitudes, once ubiquitous in particle theory, can still yield a surprising amount of in- 
formation about heavy hadron transitions. For the semileptonic decays considered in this 
paper, this stems from a two-part procedure: First, the perturbative calculation of a 
two-point function, performed in an unphysical kinematic regime where the calculation is 
reliable, is connected by crossing symmetry and a dispersion relation to the QCD form 
factors of interest. Second, some complex analysis consisting of a conformal transforma- 
tion, a multiplication by Blaschke functions, and a Taylor expansion, produces a simple 
parametrization for the form factors in the physical region. 

In the current work we have presented an improvement of the conformal transforma- 
tion that decreases the number of parameters necessary for an accurate description of the 
form factors over a given kinematic range. We have computed the explicit parametriza- 
tions for the form factors of the decay ^c^i^, and have recomputed, using the 

new conformal transformation, the parametrizations of the mesonic decay form factors for 
B D^*Uu and B iriiy. 

In particular, for a two-parameter description of A;, Aciu form factors (utilizing 
heavy quark symmetry for the normalization), we find truncation errors of 4-7% for the 
form factors associated with transversely polarized lepton pairs, and errors of 20-26% for 
the longitudinally polarized combinations. These numbers would have been much larger 
without the technical improvement to the conformal map Eq. ( 2.12| ), owing to the large 
number of Be resonance masses below the A{,Ac threshold. The truncation errors are seen 
to be quite insensitive to the exact locations of the Be poles, so our results do not depend 
on the detailed assumptions of a particular potential model calculation. 

Constraints for the decays B D^*^£u, which were already quite restrictive from 
earlier work, become more so using this technical improvement. In particular, we showed 
that one can obtain a one-parameter fit to all of these form factors except Fi good to 3- 
6%, using the normalization from heavy quark symmetry. A description of the form factor 
f+ of B —>■ Tciu good over the entire kinematic range to an absolute uncertainty of 0.37 
requires only four parameters. Relative to the expected size of f+, this represents a small 
fractional uncertainty in the pole-dominated region. In addition, if the normalization of 
at more than one kinematic point is known, the envelope of allowed parametrizations 
becomes quite restricted, even for small g^[p!6|. 
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These conclusions follow only from very general properties of QCD and ^'-matrix 
amplitudes. They do not rely on such potentially useful information as the contribu- 
tion of exited states above threshold to the dispersion relation, or the actual contribution 
(residues) of heavy resonances below threshold to singularities of the form factors. One 
might therefore expect only very weak constraints to arise; instead one is rewarded with 
strong restrictions on the shape of the form factors. In this paper, our treatment has been 
somewhat mathematical. A physical understanding of why the restrictions are so strong, 
and what role the small but apparently natural variable z plays, would be of considerable 
interest. Further detailed computations compatible with the QCD dispersion constraints 
may yield restrictions that are better still. 
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Figure Captions 

Figure 1. Important kinematic points in the z-plane for the related processes (vacuum 
Mm) and M miu using the map Eq. (|2.12|) . a) {z = —1): The branch point q'^ = 
(M + m)^, threshold for (vacuum Mm). The two sides of the branch cut > (M + m)^ 
map to the upper and lower halves of the unit circle, with q^ —>■ +00 corresponding to 
z — > +1 along the circle, b) Location of resonant poles below the (vacuum —>■ Mm) 
threshold. c,d) ^min, -Zmax defined in Eq. (|2.13|) . The segment (^min, ^max) is the kinematic 
region for M — > mii>, and is chosen in Eq. ( p.l2|) so that z = lies in this interval. 
q^ — > —00 corresponds to ^ ^ +1 along the real axis. 

Figure 2. The one-parameter description of the B D*ii? form factor g{uj)/g{l), 
plotted as a function of velocity transfer lo, for a set of values of the parameter ai. The 
lowest curve corresponds to the saturation of Eq. ( |3.3| ) by ai = —1, while the top curve 
corresponds to ai = 0.115, the largest value allowed by the the ctg— improved Bjorken 
inequality P7| , [2^ . Intermediate curves correspond to equally spaced values of ai. 
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Fig. 2: Form Factor g(w)/g(l) 




w 



